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Summary 

This report describes a methodology which predicts the 
behavior of ceramic matrix composites and has been incorpo- 
rated in the computational tool CEMCAN (CEramic Matrix 
Composite ANalyzer). The approach combines micro- 
mechanics with a unique fiber substructuring concept. In this 
new concept, the conventional unit cell (the smallest represen- 
tative volume element of the composite) of the micromechanics 
approach is modified by substructuring it into several slices and 
developing the micromechanics-based equations at the slice 
level. The methodology also takes into account nonlinear 
ceramic matrix composite (CMC) behavior due to temperature 
and the fracture initiation and progression. Important features 
of the approach and its effectiveness are described by using 
selected examples. Comparisons of predictions and limited 
experimental data are also provided. 


Introduction 

Ceramic matrix composites (CMC’s) are the subject of 
much research because new high-temperature materials are 
required for structural components and other applications in 
high-speed engines. Continuous fiber-reinforced ceramic 
matrix composites are promising candidates and efforts to 
develop them are considerable. In comparison with conven- 
tional materials, CMC’s offer several advantages: high specific 
stiffness and strength, higher toughness and nonbrittle failure 
as compared with monolithic ceramics, environmental stability, 
and wear resistance for both room- and elevated-temperature 
applications. Ceramic matrix composites are reinforced pri- 
marily to enhance toughness. A weak fiber-matrix interface 
allows the toughening mechanisms such as fiber bridging, 
crack deflection, matrix microcracking, fiber debonding, and 
fiber pullout to be combined to reduce residual stresses and 
eliminate catastrophic failures. However, the full utilization of 
these materials can only be achieved by their accurate 
micromechanical representation. 


The analysis of CMC materials requires specialized model- 
ing which considers their unique physical and mechanical 
behavior. Ceramic matrix composite materials are usually 
reinforced by high stiffness fibers , the role of which, in addition 
to providing higher stiffness, is to enhance toughness because 
the matrix material is quite brittle and fails at relatively low 
strain levels. In CMC’s the moduli of fibers and matrix are 
usually of the same order. The research in CMC analytical 
modeling under the sponsorship of HITEMP (High Tempera- 
ture Engine Materials Program) at the NASA Lewis Research 
Center in Cleveland, Ohio, has resulted in the computational 
tool CEMCAN (CEramic Matrix Composite ANalyzer) which 
predicts CMC behavior. 

The methodology incorporated in CEMCAN is based on 
micromechanics models in which customarily a representative 
volume element or unit cell is arranged in a square array pattern. 
However, the present approach employs a unique multilevel 
substructuring technique that allows the capture of greater local 
detail. This technique has four levels of substructuring: from 
laminate to ply, to subply, and then to fiber (figs. 1 and 2). The 
fiber is substructured into several slices and the micromechanics 
equations are applied at the slice level. Although the basic 
philosophy can be applied to the analysis of any continuous 
fiber-reinforced composite, the emphasis here is on the devel- 
opment of a computer code to specifically analyze and simulate 
aspects unique to ceramic matrix composites. The aspects of 
interest include varying degrees of interfacial bond around the 
fiber circumference and accounting for the fiber breaks and 
local matrix cracking that may lead to the rapid degradation of 
the interphase at higher temperatures as a result of oxidation. In 
addition, the multilevel substructuring technique used here can 
account for different fiber shapes and can integrate the effect of 
all of these aspects on composite properties and/or responses, 
in turn providing greater details about stress distribution. 

The objective of this report is to describe the methodology 
and illustrate its effectiveness with appropriate examples, which 
include computing the following: composite equivalent 
mechanical properties; microstresses in various regions due to 
varying degrees of interfacial bond; nonlinear material behav- 
ior due to temperature and stress redistribution caused by 
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Figure 1.— Integrated analysis approach embedded in 
CEMCAN computer code. 

progressive fracture and/or damage under thermal and 
mechanical loading. The predictions made, based on the present 
methodology, will also be compared with experimental results. 
Appendix A contains a list of symbols used in this report. 


Fiber Substructuring and Micromechanics 

The primary objective of composite micromechanics is to 
determine the equivalent elastic properties of a composite 
material in terms of the elastic properties of the constituent 
materials. The properties of interest are composite moduli, 
Poisson’s ratios, thermal expansion coefficients, thermal con- 
ductivities, heat capacity, and various composite strengths. 
Identifying the smallest representative unit volume element 
(RVE) is usually the starting point for any micromechanics 
approach. The micromechanics equations are derived for this 
element which is sometimes referred to as a “unit cell. It is the 
smallest region over which the stresses and strains are assumed 
to be macroscopically uniform. However, within the unit cell 
the stresses and strains are nonuniform due to the heterogeneity 
of the material. The unit cell consists of fiber, matrix, and 
possibly an interphase treated as a separate constituent. It is 
assumed that these cells are arranged in a regular square array 
pattern. Equivalent material properties for the ply are then 
derived in terms of constituent material properties based on the 
mechanics of materials approach. Other assumptions are that 
(1) the fiber and matrix are subjected to the same strain in the 
fiber direction of a unidirectional fibrous composite and (2) the 
same transverse stress is applied to both the fiber and matrix in 
the direction transverse to the fiber. These are standard assump- 
tions for a mechanics of materials approach. However, it should 
be noted that such an approach is not mathematically rigorous. 
It only ensures the force equilibrium in all directions. A 
mathematically rigorous solution that ensures the continuity of 
strains, displacements, and stresses across the fiber and matrix 
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Figure 2.— Ply/fiber substructuring concepts for 
ceramic matrix composites micromechanics. 

(a) Horizontal slicing (along 2-2 direction). 

(b) Vertical slicing (along 3-3 direction; note: 
i-1 is along fiber direction), (c) Individual slice 
showing fiber, matrix, and interphase regions. 
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boundary can be accomplished through the use of the theory of 
elasticity. However, the solution of elasticity equations for a 
composite material is quite tedious. 

The present methodology for CMC’s uses a unique fiber 
substructuring approach to derive the micromechanics equa- 
tions. The unit cell is further subdivided into several slices. The 
equations for micromechanics are derived for slices; that is, the 
slice equivalent properties are computed based on the proper- 
ties of the fiber, matrix, and interphase. The fiber substructuring 
and slice geometry are shown in figure 2. The derivation of the 
micromechanics equations and the assumptions involved with 
the approach are not given here but are documented in refer- 
ences 1 and 2. The equations, however, are presented in 
appendix B for the sake of completeness. When the equivalent 
slice properties are acquired, the equivalent properties of the 
unit cell, or RVE, are obtained by using the laminate theory in 
a manner analogous to obtaining laminate properties from ply 
properties. It should also be mentioned that 2-2, or horizontal, 
slicing is used to compute 1 - 1 and 2-2 slice properties whereas 
3-3, or vertical, slicing is used to compute slice properties in the 
3-3 direction. 

The fiber substructuring approach is general and versatile 
and can be applied to any type of continuous fiber-reinforced 
matrix composite. Also, the fiber shape or geometry does not 
pose a problem at all so that various types of fiber cross sections 
can be routinely analyzed with the present approach. Slicing 
details must be provided through user subroutines for fibers that 
do not have a circular cross section. For circular cross sections, 
the program automatically does the slicing. However, analyses 
for other fiber shapes are beyond the scope of the present report 
and therefore are not examined. By carefully controlling the 
number of slices, one can easily incorporate various degrees of 
bond around the fiber circumference and along the length, fiber 
breaks and matrix microcracking, and effects of these param- 
eters on ply thermal and mechanical properties and/or responses. 


Material Nonlinear Behavior 


At NASA Lewis over the last two decades, research to 
simulate the behavior of composite materials has culminated in 
a unified law to describe the material degradation behavior. 
This law, referred to as the multifactor interaction relationship 
(MFIR), accounts for the degradation effects of environment, 
fabrication, and load on constituent material behavior. In the 
present work, the constituent properties, namely the modulus 
and coefficient of thermal expansion, are considered a function 
of temperature only. For example the modulus is assumed to 
have the following functional relationship with temperature: 
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where E is the modulus at temperature T, E () is the reference 
modulus at the room temperature T 0 , and T f is the final 
temperature where the modulus is zero. Based on the available 
experimental data, the exponent n and the final temperature T f 
must be computed separately for fiber and matrix properties. 
Details of one such computation for the SiC/RBSN composite 
is given later in the section Results and Discussion. 


Microstresses and Stress Redistribution 
Due to Progressive Fracture 

The thermal and mechanical loads are specified or obtained 
at the laminate level. Through the successive use of lamination 
theory, the stresses and strains at the midplane of the plies and 
then at the midplane of the slices are obtained. The constituent 
microstresses are then computed from microstress equations 
that relate the equivalent slice stresses to the constituent stresses. 
These equations, which are based on the principles of 
micromechanics, are listed in appendix C. 

In the stress-strain curve (in addition to any material 
nonlinearities mentioned above), ceramic matrix composites 
also show nonlinear behavior due to matrix microcracking that 
leads to stress redistribution as the fracture initiates and propa- 
gates. These parameters can be accounted for automatically. 
The microstresses are computed in all the constituents of a slice 
first, and the stress in each constituent is examined for failures. 
In the present methodology, a maximum stress-based criterion 
is used for the initiation of damage or fracture. However, any 
other failure criterion that depends on combined stresses can be 
easily used. Accordingly, if a constituent stress has exceeded 
the corresponding strength value, the constituent is assumed to 
have failed and the corresponding modulus of that constituent 
is reduced to an almost negligible value. A slice is assumed to 
have failed in the longitudinal direction if the fiber portion of 
the slice indicated a stress higher than the fiber fracture strength 
in the longitudinal direction. Similarly, a slice is assumed to 
have failed in the transverse direction if matrix region stresses 
indicated failure in the transverse normal or shear direction, 
respectively. Once a slice fails, all the microstresses that are 
carried by that slice up to that load step are redistributed in the 
slices that have not failed. The load that was carried by the failed 
slice is then added appropriately to the laminate load and a 
laminate analysis is conducted again. The load redistribution 
due to progressive fracture is shown schematically in figure 3. 
This process is continued for a given load step until an equilib- 
rium (convergence) state is reached between the applied load 
and a damage state. Presently, the convergence is checked for 
laminate midplane strains and ply and slice strains. When all the 
strains are within 5 percent of the previous iteration, conver- 
gence is assumed to have been reached. This mechanism is 
mainly responsible for the overall nonlinear behavior in the 
laminate stress-strain curves at room temperature. A detailed 
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• Failure of slice in longitudinal 
direction governed by fiber 
failure 

• N x ' = N X + o- s i ice 1 1 x 
(slice thickness, h) 



• Slice failure in transverse 
(shear) loading governed 
by matrix/interphase failure 

• Ny' = Ny+ ff s |jce 22 x 

(slice thickness, h) 


Figure 3.— Local stress redistribution due to progressive 
fracture. Failure based on maximum stress criterion. 


description of material nonlinearities and progressive damage 
modeling can be found in reference 3. 

CEMCAN Computer Code 

The methodology just described incorporates fiber 
substructuring, material nonlinearities, progressive stress re- 
distribution effects, and micro- and macromechanics equations 
and has been programmed into a stand-alone computer code 
CEMCAN. The flowchart of the program is shown in figure 1 . 
The transition of slice properties to equivalent-unit-cell (or 
RVE) properties is achieved by using the classical laminate 
theory. If there is only one fiber through the ply thickness, the 
ply equivalent properties are identical to the properties of the 
unit cell or the RVE. If there are a number of fibers through the 
thickness of a single ply, the unit cell properties have to be 
integrated by using the classical laminate theory again to obtain 
the ply properties. Then, the laminate or composite properties 
are obtained from the single ply properties by using 
macromechanics theories. The classical laminate theory will 
not be explained here because it can be found in any composite 
mechanics textbook (ref. 2). After the laminate properties are 
computed, the next step is to obtain the laminate response due 
to externally applied loads. The various steps involved are 
depicted in figure 1. The left side of the flowchart shows the 
integration or synthesis of the properties from the slice (lowest) 
to the laminate (global) level. The right side of the flowchart 
shows the decomposition of the response from laminate to ply, 
to slices, and then to microstresses when thermal/mechanical 
loads are applied. The structure of the CEMCAN computer 
code closely resembles the previously developed composite 
mechanics codes at NASA Lewis (refs. 4 and 5). 


Currently, the code can predict the composite mechanical 
and thermal properties and can compute the stresses at ply, 
slice, or constituent levels, thus allowing the prediction of very 
detailed stress gradients at a high computational efficiency. 
Although the results are not shown herein, the thermal and 
mechanical properties, as well as the microstress predictions, 
have been verified with other micromechanics theories and 
very detailed three-dimensional finite element analyses. The 
user can specify the number of slices in a unit cell and the state 
of interfacial bonding around the fiber circumference in each 
slice; the bonding ranges from 0 (completely debonded) to 1 
(fully bonded). To obtain the mechanical properties of the 
interphase, the values provided in the resident data bank are 
multiplied by the factor (0 to ^representing interfacial bond- 
ing. Thus, if the interphase in a slice is completely debonded, 
the mechanical properties of the interphase will be reduced to 
negligible (almost zero) values. The code can also predict the 
stress-strain behavior of a laminate up to failure by considering 
the processing conditions. Two types of nonlinearities can be 
handled with CEMCAN. The first type is that arising from a 
stress redistribution in the constituents as a result of progressive 
damage of the interphase, matrix, or fiber within the unit cell. 
If a constituent fails in a slice, it unloads and its microstresses 
are redistributed to the unfailed regions. Consequently, the 
global stress-strain behavior becomes nonlinear. The second 
type of nonlinearity results from the dependence of constituent 
properties on temperature. The residual stresses resulting from 
the cooldown phase of the fabrication process can be tracked 


with the code. 

Another unique feature of the computer code is the resident 
data base of constituent properties, which is divided into three 
sections. The first, consisting of fiber properties, enables the 
user to include any number of different properties; however, 
the properties must be identified by four-character code names 
that can be used by the program to search the data base. The 
second and third sections consist of matrix and interphase 
properties. There is no specific limit on the number of different 
constituent properties one can include in the data base. The 
complete details of the code operation and the resident data base 
are discussed in the CEMCAN user’s guide (ref. 6). 


Results and Discussion 

To demonstrate the methodology and to illustrate the various 
current capabilities of the CEMCAN computer code, the com- 
posite system SiC/RBSN (silicon carbide SCS-6 fibers in 
reaction-bonded silicon nitride matrix) was chosen because of 
the vast amount of readily available in-house experimental 
data. These data are also documented in reference 7, which 
forms the basis for the computed properties that are used in the 
present study. The computation procedure used herein is 
described in the next section. The use of such a procedure is 
mandatory for any new material system in order to arrive at a set 
of constituent level (fiber, matrix, and interphase) properties 
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that will be created in the resident data base. These properties 
can then be used to predict the thermomechanical response of 
ceramic matrix composites made of any layup type. 

Prediction of Thermal and Mechanical Properties 

The chosen composite system has a very compliant and weak 
interphase. The region between the fiber and the matrix, the 
interphase bonds them and is treated as a separate constituent 
in the study. The interphase is usually a combination of an outer 
carbon layer of the fiber and any reaction zone resulting from 
the cooldown phase of the fabrication process. The interphase 
properties are not readily available and therefore must be 
inferred from other properties. Based on photomicrographs, the 
interphase thickness is estimated to be approximately 3 percent 
of the fiber diameter (0.03 x 142 pm = 4.3 pm). The nominal 
fiber volume ratio (fvr) for this composite was estimated to be 
in the range 0.3 to 0.34 in reference 7. It is well known that the 
interphase properties have very little effect on the longitudinal 
modulus of a unidirectional composite (ref. 8); therefore, we 
calculated the nominal fiber volume ratio to be 0.36 by compar- 
ing the measured and predicted values of longitudinal modulus 
for a unidirectional composite (fiber and matrix moduli are 
known). 

For the purpose of modeling, two assumptions were made: 
the formation of the interphase erodes the fiber, thereby reduc- 
ing its effective diameter, whereas the matrix material remains 


the same. Consequently, the fiber volume ratio is reduced but 
the matrix volume ratio remains unchanged. For example, in 
the present analysis, the gross fiber volume ratio was 0.36, the 
fiber diameter was 142 pm, and the interphase thickness was 
3 percent of the fiber diameter. Effectively, then, the fiber 
volume ratio was reduced to 0.3 1 8, the interphase volume ratio 
was 0.042, and the matrix volume ratio remained unchanged at 
0.64. On the other hand, if voids exist, the matrix modulus is 
appropriately reduced in proportion to the amount of voids 
internally. The user need only give the amount of voids or 
porosity. When the nominal fvr is known, the interphase 
modulus can be computed by comparing the measured and 
predicted values of the transverse modulus for a unidirectional 
composite (because the transverse composite modulus is quite 
sensitive to interfacial conditions). Alternatively, the measured 
value of the in-plane shear modulus for a unidirectional com- 
posite can also be used to compute the interphase modulus. The 
fiber and matrix properties and the computed interphase prop- 
erties are shown in table I. 

Once computed, the interphase property values can be used 
to predict the mechanical properties of laminates with general 
layups. Tables II and III show the properties predicted for three 
different laminates with fairly common layups: [0 8 ], (0 2 /90 2 ] s , 
and [±45 2 ] s . Also shown are the experimentally obtained 
properties for the same laminates. There is excellent agreement 
in the results. It should be noted that measured values are based 
on several specimens and indicate some scatter. To allow for 


TABLE I.— FIBER, MATRIX, AND CALIBRATED INTERPHASE PROPERTIES OF SiC/RBSN COMPOSITE SYSTEM 


Constituent 


Property 1 


Poisson's 

Modulus, 

Shear modulus. 

Coefficient of 

Thermal 


ratio, 

V 


E 

C 

j 

thermal expansion, 
a 

conductivity, 

K 



GPa 

Mpsi 

GPa 

Mpsi 

loVc 

ioVf 

W/m-K 

Btu/ft-hr-°F 

SiC (SCS-6) Fiber 

0.17 

390 

56.6 

117 

17 

4.1 

23 

22 

12.7 

RBSN Matrix 

.22 

110 

15.95 

45 

6.5 

22 

12 

5 

2.9 

Interphase 

.22 

3.5 

.5 

1.4 

.2 

20 

1.1 

20 

12 


TABLE II —PREDICTED (CEMCAN) AND MEASURED ELASTIC PROPERTIES OF 
SiC/RBSN f0 8 ] LAMINATE 


[Unidirectional SiC/RBSN [0 8 ] composite fiber volume ratio, 0.36 a .] 


Value 

Property 


Young’s modulus 
(longitudinal), 
E 

Young's modulus 
(transverse), 

§22 

In-plane shear 
modulus, 

^12 

Poisson's 

ratio, 

V,j 


GPa 

Mpsi 

GPa 

Mpsi 

GPa 

Mpsi 


Predicted 

189 

27.4 

68.3 

9.9 

27 

3.9 

0.21 

Measured b 

193±7 

28+1 

69±3 

10±0.4 

31±3 

4.5±0.4 

.2 


The measured value of longitudinal modulus Ejj for a unidirectional composite was 

used to compute fiber volume ratio (fvr) whereas the transverse modulus E 22 was 
used to compute the interphase modulus. 


Reference 7. 
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TABLE m — PREDICTED (CEMCAN) AND MEASURED ELASTIC PROPERTIES OF SiC/RBSN 

[02/902] and [±45 2 1 s LAMINATES 


Laminate 

Property 



Poisson's 

Young’s modulus 

Young’s modulus 

In-plane shear 


ratio. 

(longitudinal). 

(transverse). 

modulus. 


V\7 


E,i 


— 

C 

'12 



GPa 

Mpsi 

GPa 

Mpsi 

GPa 

Mpsi 

[(V90 2 ] s 

0.112 

129 

18.7 

129 

18.7 

27 

3.9 

Predicted 
Measured a 

.12 

124±6 

18±0.9 

124±6 

18±0.9 

31±2 

4.5±0.3 

[±45 2 ] s 

0.46 

78.6 

11.4 

78.6 

11.4 

58 

8.4 

Predicted 

Measured 3 

.36 

78±3 

1 1.3±0.4 

78±3 

11.3±0.4 

— 

... 


Reference 7. 


this scatter, one must systematically account for the variability 
or scatter in the properties of the constituent materials and 
perform probabilistic micro- and macromechanics analyses. 
Such procedures would generally yield the expected or mean 
values of the response (mechanical properties) as well as their 
distribution. Henceforth, one could put scatter bounds on the 
predicted values. An ongoing effort is to include a probabilistic 
description of the behavior in the methodology and to develop 
a formal computer code. 

Prediction of Stress-Strain Behavior Under Uniaxial 
Loading 

This section presents comparisons of experimental data and 
the prediction of composite behavior up to failure under uniaxial 
tensile loading at room temperature for various layups. Here 
one must account for the fabrication-induced thermal stresses 
as well as the material nonlinearities arising from the wide 
temperature range involved. To allow for the fabrication- 
induced thermal residual stresses in the composite behavior, 
one needs to specify the stress-free temperature from which the 
composite is cooled down to room temperature during the 
fabrication process. The composites are processed (i.e., 
nitridation is performed) between 1200 to 1400 °C. However, 
after performing several analyses, the stress-free temperature 
in the present simulations was assumed to be 900 °C so that the 
observed failure modes in these composites could be obtained. 
This suggests that at higher temperatures (900 to 1200 °C) some 
metallurgical mechanism is operative and prevents stress buildup 
above the assumed stress-free temperature (900 °C). Such 
observations have been made before with regard to other 
ceramic matrix composites (ref. 9). In the present simulations, 
the composites are cooled to room temperature (to 20 from 
900 °C) prior to any loading. 

The material nonlinear behavior due to temperature can be 
accounted for by using a simple functional relationship 
between the room-temperature reference properties and those 
at any temperature of the constituents. However, such an 
expression should be properly computed using available data 
for the matrix and the fiber material. The functional relationship 


used herein was given in equation (1 ). Note that the maximum 
use temperature of these composites will be about 1 200 °C and 
that very little degradation (max of 10 percent. Private commu- 
nication with Dr. RamakrishnaBhatt, Materials Division, NASA 
Lewis Research Center, Cleveland, OH, 1994.) of the constitu- 
ent modulus exists, even at the maximum use temperature as 
compared with the reference property at room temperature. To 
take that into account, the exponent for the temperature depen- 
dence of the properties given in equation ( 1 ) is 0.25 for the fiber 
and 0.1 for the matrix; the final temperature Tj is taken as 
2500 °C for the fiber and 2200 °C for the matrix material. These 
assumed values correspond to the observed degradation in the 
constituent material properties, namely the modulus and the 
thermal expansion coefficient for the use temperature range of 
these materials. Figure 4 shows the behavior of fiber and matrix 
moduli with respect to temperature according to equation (1). 

To illustrate the procedure just described, the stress-strain 
behavior up to failure for unidirectional on- and off-axis speci- 
mens and two angle-ply laminates is predicted and compared 
with experimental stress-strain curves at room temperature. 



Figure 4.— Effect of temperature on constituent (fiber 
and matrix) moduli according to equation (1), (E/Eq) = 
[{Tf-T)f(Tf-T$T- Final fiber temperature, T h 2500 °C 
for n = 0.25; final matrix temperature, T f , 2200 °C for 
n = 0.1 . 
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The details of the simulation are shown in figures 5 to 7. The 
stress-strain curves for [0 g ] and [90 g ] composites are shown in 
figure 5. Based on the observed tensile strengths of these 
composites, the matrix in situ tensile strength was calibrated to 



Figure 5. — Room-temperature stress-strain curves to 
failure of [0] 8 and [90] 8 SiC/RBSN composite. 



Figure 6. — Room-temperature stress-strain curves to 
failure of [10] 8 and [45] 8 laminates of SiC/RBSN. 



Figure 7.— Room-temperature stress-strain curves to 
failure of [±452l s SiC/RBSN laminate. 


be 1 00 MPa (14.5 ksi) by comparing the knee in the stress-strain 
curve of the [0 8 ] composite, and the average fiber bundle 
strength, to be 2 GPa (285 ksi) by comparing the ultimate 
strength of the composite. The knee represents the initiation 
and subsequent saturation of the matrix microcracking. These 
values agree well with the fiber and matrix tensile strengths 
reported in reference 7. Once calibrated, the same strength 
values were used in subsequent simulations. The stress-strain 
behavior of the [0 8 ] laminate (fig. 5) is a bilinear behavior, the 
knee of which represents matrix cracking to saturation. During 
the second part of the behavior, the load is entirely carried by 
fibers only. The behavior of the [90 g ] laminate is linear until 
failure, at which point the matrix failed in transverse tension at 
a relatively low stress level of 45 MPa (6.5 ksi). 

In contrast, the off-axis laminates [10 8 ] and [45 g J stress- 
strain behavior is largely linear up to failure as shown in 
figure 6. Such behavior is characteristic of a brittle failure. The 
final fracture in these laminates is controlled by a combination 
of matrix and interphase normal and shear failure modes. The 
room-temperature stress-strain curve to failure for a [±45-,) s 
laminate is shown in figure 7. The shear strength of the matrix 
is estimated to be half the matrix tensile strength to match the 
ultimate stress of the laminate under consideration. This lami- 
nate also fails in a combination of shear and normal failure 
modes and exhibits the characteristics of a graceful failure. 
Once again, the predicted strength of this laminate matched the 
experimentally measured value. Although the optical micro- 
graphs in reference 7 confirmed the aforementioned failure 
modes, the unloading portion of the stress-strain curve for the 
laminate was not predicted. The experimental data shown in 
these figures was taken from reference 7. 

Influence of Partial Interphase Bond 

The concept of fiber substructuring in the present method- 
ology allows one to specify a partial bond around the fiber 
circumference and then integrate its effect up to the composite 
properties and response. The variation of some mechanical and 
thermal properties as a function of the percent of the fiber 
circumference debonded is shown in figure 8. Longitudinal, or 
fiber-controlled, properties show little degradation whereas the 
transverse properties reveal greater degradation as a function of 
debonding. Figure 9 is a comparison of composite properties 
normalized with respect to the value of the property for a 
composite with a strong interphase for a varying degree of 
interface bonding through the thickness. The layup is [0 5 ]and 
the fiber volume ratio of the composite is 30 percent. The figure 
considers three levels of interfacial damage that are assumed to 
have taken place as a consequence of oxidation damage from 
exposure to high temperatures. The first case is interfacial 
damage only in the outer plies whereas the second case is a 
linearly varying degree of interfacial damage with maximum 
damage in the outer plies. In the third case, all the plies have the 
same degree of interfacial damage. The figure reveals that the 
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Percent circumference debonded 

Figure 8. — Variation of thermal and mechanical 
properties (normalized with respect to property 
with strong interphase) due to partial bonding. 


Damage 


High *** 



£ C 1 1 • "cl 1 E c22. E c33 G c12. G c23 *VI2 a c22 

Figure 9. — Comparison of composite properties (normalized with 
respect to property with undamaged interphase) due to interfacial 
damage through the thickness; fiber volume ratio, 0.3. 


damage in the interphase region leads to a varying degree of 
bonding and more severely effects the transverse stiffness 
behavior. There is a substantial reduction in the composite 
transverse normal and in-plane shear moduli; the longitudinal 
behavior is relatively unaffected. 


Conclusions 

A unique and novel fiber substructuring technique has been 
incorporated in the computer code CEMCAN to predict the 
behavior of ceramic matrix composites. Its features are illus- 
trated with selected examples. Experimental validation is pro- 
vided where the data are available. This concept enables one to 
study in great local detail ceramic matrix composite behavior 
pertaining to interfacial bonding of varying degrees. Based on 
the results shown for the SiC/RBSN composite system (silicon 


carbide SCS-6 fibers in reaction-bonded silicon nitride matrix), 
the following specific conclusions can be drawn: 

1 . Fiber substructuring captures and represents greater local 
detail than other unit-cell-based micromechanics theories. 

2. CEMCAN predictions have been verified successfully 
for the SiC/RBSN composite system for which experimentally 
obtained data are available. 

3. The material nonlinear behavior model and distribution 
of progressive damage with stress have enabled the prediction 
of stress-strain behavior up to failure for different laminates. 
The agreement with experimentally observed behavior is 
excellent for the SiC/RBSN ceramic composite system. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, April 26, 1996 
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Appendix A 
Symbols 


c 

heat capacity 

p density 

E 

Young’s modulus 

o stress 

G 

shear modulus 

Subscripts: 

K 

thermal conductivity 


k 

volume ratio 

c composite or laminate 

T 

temperature 

/ fiber 

a 

coefficient of thermal expansion 

i interphase 

£ 

strain 

i slice 

V 

Poisson’s ratio 

m matrix 



Appendix B 

Micromechanics Equations 


The micromechanics equations describing the relation 
between the constituent and the equivalent ply properties are 
derived for the unit cell using the mechanics of materials 
approach. The complete details can be found in reference 1 . 
Only a brief description of the equations is given here for the 
sake of completeness. The modeling details for the substructured 
unit cell and the slice are shown in figure 2. 


Equivalent Slice Properties 


To derive the equivalent slice properties, consider a slice 
taken from the unit cell. Let dj , d m , and d i be the fiber, matrix, 
and interphase widths and h and s be the height of the slice and 
the total width of the slice, respectively. Because it is customary 
to express the equivalent properties in terms of the properties of 
each constituent and its respective volume ratio, let kpk m , and 
k i be the fiber, matrix, and interphase volume ratios, respec- 
tively. Then, by definition. 



2 d m . _ 2 d, 

S 5 


(Bl) 


Also, by definition 

\ jn = (B6) 

£ ;ll 

where j can represent any of the subscripts m, or /. 

The combination of equations (Bl ), (B5), and (B6) leads to 
the following equation for the major Poisson’s ratio: 

v m ~ k f V f\2 + k m v m\2 + *j v i'12 ^ B7 ) 

The modulus in the transverse 2-2 direction can be derived by 
assuming that the same transverse stress a f22 * s a PP^ e ^ to the 
fiber, matrix, and interphase: 


<*£22 ~ °m22 ~ °i22 “ °/22 (®8) 

With the aid of equations (5) and (8) and noting that one can 
derive the following equation for the transverse modulus in the 
2-2 direction 


Now that the preliminary parameters have been established, the 
mechanics of materials approach is applied to generate the 
composite mechanical and thermal properties. 



(B9) 


Mechanical Properties 

By applying the force equilibrium in the longitudinal direc- 
tion, the following equation can be written: 

o ni hs = (2o m i]dm +2o in d i +Ofnd/)h (B2) 

Also, based on the assumptions for the mechanics of materials 
approach, the following equation can be written as 

e m ~ e mll = e il l = e /i l 


fc *22 


E m22 E i22 E j22 

k m E i22 E f22 + k i E f22 E m22 + k f E m22 E i22 


(BIO) 


The in-plane shear modulus in the 1-2 direction is determined 
by assuming that the shearing stresses on the fiber and the 
matrix are the same. In this respect, the derivation of G m is 
similar to that for 2 * The equation for G n2 


G m ~ ' 


G m\2 G i\2 G f\2 


k m G i\2 G f\2 + k i G f\2 G m\2+ k fG m \2 G i\2 


(Bll) 


With the aid of equations (B 1 ) to (B3), the longitudinal modu- 
lus E n j can be written as 


The remaining mechanical properties can be derived in a 
similar fashion and are given by 


E n \ - k m E m\\ ~ k i E n\ (B4) 

The Poisson’s ratio v in the 1-2 direction can be derived by 
considering the strain in the transverse direction 2-2 as the sum 
of the individual strains in the constituents. Accordingly, 

E ?22 s = ^ z m22 d m + ^ z i22 d i + z f22 d f (B5) 


E m - 


E m 33 E i 33 E f 33 


k m E i 33 E /33 + k i E f 33 E m 33 + k f E m 33 E i 33 


(B12) 


G m - 


G m\ 3 G i\ 3 G f \3 


k m G iU G f \3 + k i G f\ 3 G m \3 + k f G m\ 3 G i \3 


(B13) 
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G m - - 


G m23 G i23 G f23 


(B14) 


k m G i23 G f23 + k i G f23 G m23 + k f G m23 G i23 

v m = */ v /13 + *m v ml3 + *i v i“13 ( B15 ) 

V 02 = k f V f32 + k m v m32 + */ v 32 ( B16 ) 


The equivalent thermal conductivities in the 2-2 direction can 
be derived by noting that the total heat flow rate across each 
constituent in the 2-2 direction should be the same. Thus, 
figure 10 shows that T 0 through T 5 are the temperatures at 
various cross sections along the slice. Equation (B21) can be 
rearranged as 


by the reciprocity relations 


*£23 


'*21 




_ v 02^122 


v 11 2 ^122 

E n\ 

E n\ 


CB19) 


given 

Q = K m (Lh) 

1 1 

Jr* 

^ i 

= K m 22 (Lh) 

1 1 

J? 

1 J 

(B17) 

= K,22 (Lh) 

1 1 

2^ 
^ i 

i i 

= Kf 22 (Lh) 

1 1 

1 + 
e^t 

(B18) 

= K i22 (Lh) 

f 

1 

i 

= K m22 (Lh) 

i 

h?: 

1 

h^| 


4 J 

d m J 


(B22) 


Thermal Properties 

The equivalent thermal conductivity in the 1 - 1 direction can 
be derived by equating the total heat flow rate across the cross 
section to the sum of the individual rates in each constituent. 
The heat flow rate can be calculated using the Fourier’s law of 
heat conduction. By following the previous steps, the heat- 
flow-rate equilibrium equation can be written as 

„ u dT , , i dT , . dT v , . dT 

K n k f\ 1 dfh — + K m\\ d m h — + Ki\ \ d i k ~T 

dL dL dL dL 

(B20) 

where K is the thermal conductivity and ( dTidL ) is the thermal 
gradient across a typical length of the cell L in the 1-1 direction 
(fig. 1 0). Simplification of equation (B20) leads to the following: 


A^n - kfKf n + k m K m n + 


ll 


(B21) 



Figure 1 0. — Modeling details of slice for deriving 
equivalent slice thermal conductivities. 


K m ~ K m22 


t,-t a 

t 5 -t 0 


s 

dm 


t 5 -t 4 


(r 5 - t 4 ) + (t 4 - r 3 ) + (r 3 - t 2 ) + (r 2 - r, ) + (r, - 7b) 

(B23) 


By observing that (7 4 - T 3 ) = (T 2 - T } ) and that (T 5 ~ T 4 ) 
= (Tj - T 0 ), equation (B23) can be further simplified to 


a: 


*22 


2 + 2 


V T 5 - T A ) 


h-T 2 

V T 5 ~ T 4 ) 


\ K m22 ' 


(B24) 


With the aid of equations (Bl) and (B22), equation (B24) can 
be written as 


V 122 


v m22 


which leads to 




k , Mm2 2 k f K m22 

m *i22 K f22 


KjrinKj 22^/22 


Kill K fll + kjKf 22 K m2 2 + kfK m22 K i22 


(B25) 


The thermal conductivity along the 3-3 direction can be derived 
in a similar manner and is given by 
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K m - 


K m33 K i33 K f33 

^m^/33^/33 + ^£^/33^m33 + */*m33^/33 


where AT is the temperature difference. Simplification of 
(B26) equation (B29) leads to the following expression for the longi- 
tudinal thermal expansion coefficient: 


The equivalent heat capacity C t can be derived by equating the 
sum of the heat contained in the individual regions to the heat 
contained in the slice based on the equivalent heat capacity: 

Q — ^"Pm^vn ^/P/Q ^/P fCf (B27) 

' ~ P t 


k m a mU E m\\ + *i a ill £ ill + k f a f\l e fU /D , m 

a ni “ / (B30) 

L n\ 

The transverse thermal expansion coefficient can be derived by 
equating the total strain in the 2-2 direction to the sum of the 
individual strains in each region: 


where p f is the equivalent density of the slice and p m , p 1? and 
p , are respectively the densities of the matrix, interphase, and 
fiber. The equivalent density can be expressed in terms of the 
individual densities by the simple rule of mixtures: 


~ ^m a m22^ + ^i a i22 + <*/ a f22^ (** 31 ) 

With the aid of equations (B 1 ) and (B3 1 ), a can be written as 

a (22 =k f a f22 + k m a m22 + k, a i22 (B32) 


p f =kfpf+ k m p m + kfii (B28) 

The longitudinal thermal expansion coefficient a m can be 
derived by noting that the sum of the forces in the longitudinal 
direction should be zero: 

{ 2d m E m\] + 2d i E iU + d f E fu) a (\] AT 

= -( 2d m a mU E mU + 2d i a iU E i\\ + d f a f\] E fll) AT ~ G 

(B29) 


The expression for oc^ 33 is similar to equation (B32) and is 
given by 

a m = k f a /33 + * m a m33 + *,a l33 (B33) 

This set of equations defines the equivalent properties for a 
typical slice with three regions. The procedure to obtain the 
equivalent properties of the representative volume element 
(unit cell) is analogous to that of obtaining the properties of a 
typical laminate from the ply properties. Here, each slice is 
treated as a subply within the unit cell. 
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Appendix C 

Equations for Microstresses in Slice Regions 


Microstresses Due To Applied Mechanical Loads 

Following are the microstress equations due to applied 
mechanical and thermal loading. 


O/n =®m + 0 ( 11 + °m 


Matrix Microstresses Due To Applied Slice Normal Stresses 


„(2) -rr 
CT i22 — ^ (22 


a (D _ E m\\ a 

L tU J 


o = 0 


„(2) _ ( v v \ a n2 F 

a mll - _ l V #21 “ V m2lj— E n 


a 54=0 


°mll - ( v /3 1 " v m3l)'7T^' £ mll 


The case for 3-3 is the same; just replace 2-2 with 3-3. 

Fiber Microstresses Due To Applied Slice Normal 
Stresses 


rr -rr (1 > 4- rr (2 > + rr (3 > 

°m\] - CT mll +G mll + G mll 


O o) .r/ 11 L 

c f it - a ^ii 

Zen 


n (2) — n 

°m 22 ~°n 2 


Poisson’s effect is computed from equilibrium considerations: 




+ k a ^ + k a ^ — 0 
V a /ll +/c i a /ll ~ U 


< 22=0 

The case for 3-3 is the same; just replace 2-2 with 3-3. 

Interphase Microstresses Due To Applied Slice Normal 
Stresses 


o 0) _ E iW L 
CT ni ~ 77 °dl 


CT iH - - ( v *21 _ V (2l)‘7^' E ill 
n (22 


Microstresses Due To Applied Slice Shear Stresses 


°/12 ~ °/12 ~ CT ml2 “ a (l2 


ct /13 - °il3 ~ °ml3 ~ °£13 


ct /23 ~ °/23 - a m23 ~ ° <23 


Microstrains 


°ni =-( v f3i- v «3i)^ £2i E,i I 


_ _ °/12 . _ _ <*,12 . _ _ <* ml 2 

e /] 2 -t; - £, 12-7; — • £,7,12-7; — 
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Microstresses Due To Applied Thermal 
Loads 

Matrix Microstresses 

a mll “( a m ~ a m\\) ATE m\\ 
a m22 = ( a ply 22 “ a ^22 ) ATE m22 

Interphase Microstresses 

°m = ( a ni ~ a n\) ATE i\\ 


a /22 = ( a ply 22 ” a f.22) ATE i22 

Fiber Microstresses 

a /n = ( a *n - a f\\) ATE f\\ 

a /22 = ( a p]y 22 ~ a n2) ATE f22 

Note: Restrained stresses/strains will be generated for the unit 
cell through the lamination theory application to the 
substructured fiber. 
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